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ABSTRACT 

An investigation is made into the accuracy of the description of the 
actual mixing-ratio profile by a power-law approximation. The approxima- 
tion is obtained by a leas t- squares technique in order to derive the best- 
fit exponent X of the pressure profile. The X -values were computed 
for soundings at 12 geographically diverse stations over the United States 
during the period 16-23 March, 1971. The X -values at each station 
were found to undergo a time-variation based mainly on the synoptic-scale 
variations in the mixing-ratio profiles. A stepwise multiple regression 
procedure involving up to four variables from the temperature-humidity 
soundings was utilized in order to "predict” the value of X from gross- 
parameters of the soundings by station. Tests were performed at four 
stations in comparing the observed precipitable water vapor W Q with that 
computed from the power-profile W cs which depended upon X computed from 
the sounding. The values of W cs gave a high correlation, R(W CS , W Q ) 0.97, 
for each of the four stations. However, the values of W cm using X from 
the predictive multiple regression equation gave somewhat larger standard 
errors at each of the four stations. The W cm -values when compared to W Q 
indicated some measure of predictive skill. 
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I. INTRODUCTION 



Reitan (1963) established a correlation factor of 0.98 between the 
logarithm of the monthly precipitable water vapor and monthly dewpoints 
based on 540 observations at numerous locations. The regression formula 
used by Reitan was 



In W = A + B t 



(i) 



where A = -0.981 and B = 0.0341 are the best fit coefficients resulting 
from the regression procedure. 

Smith (1966) has shown that an equation of form (1) may be derived 
from the well known empirical equation of Tetens (1930) 

e = E q x 10 (a t d " B ) / ( t d + 7 ) (2) 

provided that the mixing-ratio profile may be parameterized. The constants 
are E 0 = 6108 dynes cm“^, a = 7.5, B = 238.1 F, and Y = 395.1 F. The 
required parameterization of the mixing-ratio can be accomplished with 



reasonable accuracy by the power law approximation 
W = W Q (p/PoD 

By combining equations (2) and (3) with 

w o = £ ( e 0 /P 0 > 

Smith was able to show that 

In W = In f E Q /g + 2.3026 (a t^ - B)/440 



(3) 



(4) 



(5) 



- . . . -ln( A + 1) 

Since climatological values of the exponent A are variable with latitude. 
Smith concluded that the coefficient A of (1) is likewise dependent upon 
latitude in Reitan' s formulation. 
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TABLE I 



\ -values from Smith (1966) 



LATITUDINAL 










ANNUAL 


ZONE -°N 


WINTER 


SPRING 


SUMMER 


FALL 


AVERAGE 


0-10 


3.37 


2.85 


2.80 


2.64 


2.91 


10 - 20 


2.99 


3.02 


2.70 


2.93 


2.91 


20 - 30 


3.60 


3.00 


2.98 


2.93 


3.21 


30 - 40 


3.04 


3.11 


2.92 


2.94 


3.00 


40 - 50 


2.70 


2.95 


2.77 


2.71 


2.78 


50 - 60 


2.52 


3.07 


2.67 


2.93 


2.79 


60 - 70 


1.76 


2.69 


2.61 


2.61 


2.41 


70 - 80 


1.60 


1.67 


2.24 


2.63 


2.03 


80 - 90 


1.11 


1.44 


1.94 


2.02 


1.62 


NORTHERN 

HEMISPHERE 

AVERAGE 


2.52 


2.64 


2.62 


2.70 


2.61 
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Bolsenga has attempted to adapt Reitan’s formulation (1) to daily 
and even to hourly observations and finds that the explained variance of 
total water vapor falls off to 72% and 64 % in the respective cases. This 
could have been anticipated since total water vapor depends in each case 
upon the entire mixing-ratio profile, which depends less upon surface 
dewpoint when a decreased time interval of sampling is employed. 

The equivalence of equations (1) and (2) to the X -profile (3) holds 
best, according to Smith, when applied to climatological data. Climato- 
logical values of X were calculated by Smith (1966), from London’s global 
mean radiosonde data (1957), with resultant X -values presented as a 
function of latitude and season as shown in Table I. 

The concept of an accurate mixing-ratio moisture profile from sound- 
ings composed of few or widely spaced levels is a very valuable tool. 

It is especially important in the field of radiation flux measurement 
where values for the mixing-ratio (moisture profile description) are 
needed at all levels to the top of the atmosphere but are frequently not 
available from sounding levels much above 500 mb. Recently Plante (1973) 
and Martin (1973) have extended the use of Smith’s power-law profile to 
one of the form 

w = w 500(p/ 500 )* (6) 

with the reference or key level at 500 mb rather than a surface reference 
level which has a diurnally variable mixing ratio. Martin (1973) using 
the Fleet Numerical Weather Central, Monterey, California, computer print- 
out of mixing ratio at six pressure levels in the vertical at FNWC grid 
points found a correlation coefficient of .998 or higher between observed 
mixing-ratio and a mixing-ratio computed by Eq. (6), based on X -values 
calculated at each grid point for the 75 grid points considered. Martin 
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concluded that the power-fit as described by Eq. (6) is accurate especially 
with soundings having widely separated reporting levels in the vertical. 

In this thesis the following questions will be investigated: for a 

given station for a given period of time, can a power-fit of the nature 
of (6) be obtained which is representative of the actual moisture profile 
at 50 mb-level increments, and what degree of accuracy can be expected? 
Further, the conditions which would tend to invalidate the accuracy of 
the method will be delineated. Two methods of verification of the results 
will be considered: 

(1) the explained variance of mixing-ratio resulting from the use of 
power-law profile, and 

(2) the accuracy for the calculated precipitable water obtained from 
the profile-parameterization compared to that observed. 
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II. DATA PROCESSING 



A. THE ORIGINAL DATA 

For the water vapor analysis, twelve stations were selected on the 
basis of their locations and the availability of data. The time chosen 
was an eight day period from 16 March to 24 March, 1971. The station 
locations, elevations, and mean surface pressures are shown in Table II. 

As can be seen from Figure 1, the stations selected lie in a band across 
the southern, and in a second band across the northern continental United 
States. The data was obtained from NWSED, Ashville, North Carolina. All 
data levels were selected at levels corresponding to integral multiples 
of 50 mb (up to 400 mb), in addition to the surface level. 

The following analysis was used as the basis of solving for A , 
starting from Eq, (6): 

log w(e)_ = X log _p (7) 

W500 oDU 

Here p ranged from the surface value to p = 400 mb, and corresponding 
values w(p) were computed at each level p using sounding- level data. If 
the left side of Eq. (7) is denoted as Y, then this equation becomes 

Y = AX (8) 

where X = log (p/500). In the average sounding which extends to 
p = 1000 mb, there would be 14 simultaneous values (Y^ , X-^) so that only 
a best-fit solution for A is feasible. 

This amounts to solving the linear regression equation, Eq. (8), for 
the coefficient A of the independent variable, with the condition of the 
zero- intercept forced on the constant term. In the program library of the 
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Figure 1: Locations of selected stations. 



W. R. Church Computer Center, the most convenient program for performing 
this objective was BIMED 02R (Dixon, 1966). While this program is primarily 
used for stepwise regression in the case of multivariate regression analysis, 
it also has the desirable option of making the required transformations of 
variables first to ratios, and then to their logs as in Eq. (7). This 
requires that the sounding data has been prepared in properly formatted 
form on IBM cards. 

B. CALCULATION OF MIXING-RATIOS AND OF \ 

Having decided on the method of handling the regression, the next 
computation was that of mixing-ratio to be specified by the power-profile 
fit, Eq. (6). 

The data as it was received (by courtesy of Mr. Russell Schwanz) , was 
card-punched in a form containing the pressure, temperature (C) , and 
relative humidity at 50 mb-intervals from the surface to the top of the 
moisture sounding (which was generally at 400 mb) . The next step then 
was to compute the mixing-ratio in its approximate form 

w = .62197 (RH) e s (T) (9) 

P 

where e s is the saturation vapor pressure computed over a plane water 
surface for T 2 0 C, and e g is computed with respect to equilibrium 
over a pure, plane ice surface at temperatures T < 0 C. In doing this at 
each sounding level, the two versions of the Goff-Gratch equations for e sw 
and e s i, respectively, were used in the forms indicated below as Eq. (10) 
and Eq. (11) . 

logio e S w = -7.90298 (T s -1) + 5.02808 logio (T s )- ... (10) 

TT T~ 

... -1.3816 x 10-7 (loll • 334 (1-T/Ts) -1) + ... 

... +8.1328 x (10 - 3.49149(T S /T-1) -i) + logio ews 
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lo S10 e si = -9 . 09718(T o /T-l) -3 . 56654 log 10 (T o /T) + ... 



(ID 



... +0.876793 (1-T/T Q ) + log^Q e io 
where T s = 373.16 K and T q = 273.16 K. 

In (10) e sw is the saturation vapor-pressure of pure water at the steam- 
point temperature and e s ± is that of pure water at the frost-point 
temperature (T < 0 C) . Then w(p) is evaluated from either equation (10) 
or (11) utilizing Eq. (9). Once the values for w were obtained at each 
level in each sounding the best-fit value for A was obtained in the 
manner suggested in connection with Eq. (8) and Section III (A) . Thus a 
unique value of A resulted for each sounding. The A -values for each 
station were then grouped with other synoptic sounding data for each 
individual sounding in order to account for the reduced variance of A 
for the sounding- sample for the particular station. In this connection, 
the stepwise regression program BIMED 02R was used to determine the 
explained variance of A at each station. Comparisons were made of 
computed w-values resulting from Eq. (6) versus observed w-values over 
the sounding station. 

C. CALCULATION OF PRECIPITABLE WATER VAPOR 

Precipitable water was obtained by integration for both the observed 
soundings and for the corresponding A -profile cases and a comparison 
was made for several stations. The formulas used in the calculations are 
shown as follows. 

a. The observed value of precipitable water vapor is 

W = / w dp = 1 r V wi + W 2 A P + W 2 + w-j A P +. . . ( 12 ) 

J g 980 LA 2 1.5 2 2.5 

P 1 
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. . . + w 12 + w 13 A p + w 13 + A p 

2 12.5 2 13.5 



W 



b. The model value of precipitable water vapor is 
Ps 

/ w 500 / -2-\ I dp = w 500 (5QQ) 

\ 500/ g g (A + 1) 

Ptop 




Ps 

400 



= M500 (500) \(M A +1 -/400Y +1 
980 ( A + 1 ) I y500/ \50o/ 



( 13 ) 



The comparison of the sample sets of the two precipitable water values 
will be used as a second measure to determine the accuracy of the \ - 
profile method. 
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III. STATISTICAL CONCEPTS 



A. THE LEAST SQUARES FIT 

In this thesis, the least- squares fit was employed in at least two 
contexts. First, it was employed in solving Eq. (8) for A at each 
station and sounding time. Here the variables Y and X are given point- 
wise by 



X. = log P and Y. = log w. 

500 1 



(14), (15) 



w 500 



with i representing a level numbered from one at the surface, to two at 
the next higher level (where P first becomes an integral multiple of 
50 mb), etc. Standard texts in statistics (e. g. Brownlee, 1960) show 
that the parameter A of the power-profile as in Eq. (8) can be deduced 
from the best-fit condition 



A 




(16) 



If the zero-intercept . had not been required in (16), A would have been 
estimated from 



A = 1 E(Xj - X)(Yj - Y) 

N 

and the constant term of the best-fit simple regression would 



(17) 

then be 



A = 




(18) 



where N is the number of sounding- level observations entering into the 
power-profile (normally here N = 14). 
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The combination of conditions (17) and (18) gives the standard best- 
fit solution for the simple independent variable case 

Y = A + A X (19) 

Equation (16) gives the corresponding result when the zero-intercept 
(Y-^ = = 0) is required. In either case A. is related to the simple 

correlation coefficient R]_2 and the variables Y and X through the standard 
formula 



R 12 = 




( 20 ) 



Hence a value of Ryjx c l° se to unity means an effective least-squares fit 
by the power-profile and simultaneously a reliable value of X .In 
Eq. (20), the statistical parameters O ^ and o ^ are defined as 






2 2 2 
Yi + Y 2 + . . - Y N 



b. 



*x = 



( 



X^ + x\ 



N 



+ 



)' 

)' 



( 21 ) 



Note that ° y and o ^ of (21) are applicable to the zero -intercept case 
and apply to the more general case of the standard regression formula 

Yi - Y = R 12 _fY (X 

a X \ 

when Y^ is replaced by Y- - Y and X^ is replaced by X^ - X. 

It was possible to apply BIMED 02R regression procedure to the 
14-values of X-^ and Y^ per sounding (that is, with each sounding listed 
on 14 cards) and to generate not only the \ -value for that sounding, 
but also the simple correlation Rq2> w hich is a measure of the ’’goodness 
of fit" of the power-profile for the particular sounding case. Such 






( 22 ) 
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results for \ and R^2 h ave been tabulated under each station and for 
each observation time (Tables III through XIV). 

B. THE STEPWISE REGRESSION PROCEDURE 

If there are N samples of the statistically related variables 

. . . ,X^) i = 1,...,N, it is possible to apply the least- 

squares concept to minimize the least-squares form 

k 

^ (^li ” Ai - “ • • • "^k^ki ) minimum (23) 

i = l 

The procedure of selecting the least-squares fit of coefficients to 
minimize Eq. (23) is called multiple linear regression. In this process, 
there is no assurance that the variables X 2 ,....,X^ used to explain the 
variation of Y-^ are the most effective, and, in fact, some of the X^ may 
not contribute to the explained variance in any significant manner. 

BIMED 02R is a statistical procedure for developing a multiple 
regression equation in which the variable X^ added at each step accomplishes 
the greatest reduction in the unexplained sum-of - squares , X (Y.-Y) 2 . If 
k-1 variables are to be added in a stepwise manner the first one added 
is that which has the largest simple correlation coefficient R^ of all 
k-1 possible comparisons. As a result, the first estimate to the desired 
multiple regression equation may be obtained in the form: 

Yf = an + a 12 X 2 (24) 

with a-Q, a -^2 selected as in (17), (18). A property of this first-step 
selection is that the sum-of-squares explained by Eq. (24) is given by 

(Y11-Yi5+...+(Yin -Yi) 2 J (25) 

A result which is proved in standard statistics texts (e.g. Crow et al, 
1960) . 



(S.S. Expl) i = R 12 
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The BIMED 02R program then forms the data-sample residuals . After one 



step, all values have the form 




(26) 



( 1 ) 



The program then examines all of the simple correlations of Yl on 
Xk (the remaining independent variables) and selects that which has the 



highest simple correlation coefficient among R.(Y^ jX^), i = 3,...,k. 

Since this simple correlation is based upon the residual after the effect 
of variable X 2 has been removed (see Eq. (26) ), it is typically called 
the first order partial correla tion with the notation ^ 13 . 2 ? the under- 
standing being that the variable having this highest correlation at 

step 2 and thus being selected is called X 3 . Then the regression equation 

( 1 ) 

relating Yj_ to X 3 is obtained as described in connection with Eqs. (21), 
(22), (23) and assumes the form 



As a result of the second-step simple regression, Eq. (27), the sum-of- 
squares of Y-^ (not Y^ ^) is now statistically explained in the amount 



The statistical parameter R-^ 23 i s the so called multiple correlation at 
step 2, that is to say of variableY-^ depending jointly on the variations 
of both X 2 and X^. These independent variables were selected according 
to the criteria of having the largest simple correlation, and the largest 
first-order partial correlations respectively. 

The multiple correlation coefficient satisfies the important chain- 
rule of multiplication 



(i) 




(27) 




(28) 




( 29 ) 
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and after k -1 steps, this rule assumes the form 



1 ' R 1 . 23 * * • k “ ( L " r 12^ ' R L3.2* • -k^ • • • ^ " R Lk. 23" • 'k-0 ( 30 ) 

2 

Just as 23 describes the fractional explained sum-of -squares of Y in 
2 

Eq. (28), so 1-R^ 23 describes the fractional unexplained sum-of - squares 

about the regression surface. This parameter has decreased relative to 

2 2 

that of step one, provided R^ 23 > ^ 12 ‘ The standard error of estimate 

is related to the multiple correlation coefficient by 

2 2 2 x 

Syj x " ° y (^”^ 12 ) (31) 

The parameter standard error, Sy|x> will be listed in Tables III,..., XIV. 

In (31), the subscript YIX denotes the best fit of Y on all X^. 

The most attractive feature of the stepwise regression is embodied in 
the chain-rule of Eqs. (29), (30), with each step making the selection 

of the highest partial correlation coefficient for use in Eq. (30). The 
result is that when the unexplained fractional sum-of -squares does not 
become significantly smaller, the process of adding variables to the right- 
hand side of 

Y l = A x + A 2 X 2 + A 3 X 3 + ... A k X k (32) 

should be terminated at the (k-l)th step. 

In Section IV, application of BIMED 02R will be used on the 16-sample 
values of X derived at each individual station. This is done since it 
turns out that a time variation in X occurs with the synoptic situation 
at individual stations. In some cases this variation, as indicated by 
the standard deviation a ^ , is so small as to be limited by the accuracy 
of the data. In other cases, a rather large variation of X occurs and 
it was desired to consider how the time-varying structure of the soundings 
contribute to the X -variation. For this further study each X -value 
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was added on an IBM card containing certain relevent parameters of the 
sounding from which A was derived. For all soundings and stations a 
set of six independent parameters (for example T(850) and w(850) and 
others to be described in Section IV) were included along with the A - 

value for the identical sounding. 

Thus BIMED 02R was also used in the important stepwise sense in 
describing the A -variation in the sounding structure over the set of 
the 16 soundings per station. The detailed discussion of the variables 
used in the regression for A is given in Section IV (B, C) , together 
with the discussion of the results. 

C. ESTIMATES OF STATISTICAL SIGNIFICANCE 

1 . Of A -Profiles at Sounding Stations 

Here the correlation coefficient R-^ of the regression equation 
(7) is the basis of the estimation of the accuracy of the A -profile. 
The significance of the test is based upon Eq. (25) which gives the 
result 

2 

Fractional Explained Variance = R-^ (33) 

whereas the fraction residual (or unexplained) variance is given by 

2 

Fractional Unexplained Variance = l-R^ (34) 

The significance test is based upon the F-statistic (which is listed in 

the BIMED 02R output) 2 

£ 

F(l, 13) = 12 (35) 

(1-R2 )/ 13 
12 

The critical value of F(l,13) at the 99 % confidence level is listed in 
statistical tables as 
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F. 0L (1> 13) = 9.07 



This 99% confidence level is satisfied whenever R^ of Eq. (34) exceeds 
Rf 2 (critical ) ^ .641 

In all but one of the X -profiles, R^ exceeded this value, and in this 

exception R^ was 0.614, corresponding to a confidence level of 97.5%. 

2 

In summary, the statistic R^ listed for each sounding at the 12 
stations, is the explained fraction of the variance of Y in that sounding 
(see Tables III, . . . , XIV) . 

Finally, BIMED 02R lists also the standard error of the coefficient \ 

Standard Error of A = A *1 l-R-^ ; N = 14 (36) 

' R 2 l 2 (K-l) 

This statistic appears in column 3 of Tables III,... XIV, and describes 
the sharpness of the \ -fit in Eq. (7). 

2 . Significance Tests of the Temporal Variation of X 



Here the fractional explained variance may be computed from 



Eq. (28) after the kth entry of an independent variable X drawn from the 

k 

synoptic sounding variables has been included in 



A = A q '+ A-^X^ + . • • + A^X^ (37) 

The fractional explained variance is given by R^ in accordance with 
(28), whereas the fractional residual variance is 1-R^. In general when 
(A) is sizeable, it is possible to generate a multiple regression equation 
of form (37), which significantly reduces the residual variance. In 
general the criterion used for entering additional variables for explaining 
the X -variance was based on a stepwise selection which added at least 
2% to the total explained fractional variance. This strategy was used 
because the F-statistic for this case is 
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F k (1, N-k-1) = 



(38) 



2 

(1 - R k )/(N - k - 1) 

and here the critical value of is difficult to surpass with N 
observations. The added fractional explained variance at step k 
identical to the numerator of Eq. (38), and is required to be at 
27o at each step. 



= 16 
is 

least 
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IV. RESULTS 



A. DISCUSSION OF X -PROFILES 

This particular phase of the research began with computer cards 
containing the following individual sounding data: pressure (mb), tempera 

ture (K) , and mixing ratio (gm/kg) . There were 16 soundings of up to 14 
levels at each of the 12 stations (Table II). 

The BIMED 02R program was utilized in a least-squares context to 
solve Eqs. 7 and 16 for a best fit X for each sounding at each station. 
The " zero-intercept" was chosen in the correlation so that the computer 
solution to Eq. (7) was of the form of Eq. (8) instead of Eq. (19). A 
discussion of the statistical methods and parameters is given in 
Section III (A). Table III lists the best-fit X -value for each avail- 
able sounding at Waycross, Ga . , over the test period. 

It should be noted that the least-square procedure for determining 
from 

Y. = A X. + e . (39) 

i'ii v 

2 

minimizes the sum-of -squares of the residuals Hence the estimator 

Y may be defined as 

Y = A log P/500 = log w c (P) (40) 

w 500 

where the observed w^qq appears in the denominator of Eq. (40) by the zero 
intercept concept. On the other hand, the observed Y is defined as 

Y = log w/w^qq so that the residual becomes 

£ i = Y - Y = log w c (P) (41) 

w(P) 
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In Table III, successive values of X obtained by X -profile fit 
[Eq. (39)] have been listed along with the standard error of each X 
determination. The small standard error of X indicates the close fit 
of the X values by Eq. (16). 

2 

In addition, the fractional explained variance R-^ for each sounding 

is listed alongside R-^ under column 4. Finally, the standard error of 

estimate is listed in the final tabular column as defined by Eq. (31). 

The value of Sy|x i- s related to the standard deviation cTy through 

2 

^YIX = (1 ~ ^12^ a X (42) 

These ratios of Eq. (42) obtained by applying the values of S YI X and 
from columns 5 and 6 of Tables III through XIV are generally of the 
order of 0.3 or lower. There are exceptions, but these can be identified 
by the values of R^ (generally R^ 9 > .95) for the respective soundings. 

In Table III, the standard error may be written in terms of the mean 
residuals - squared as 



r 


N 


/ \ 0 




1 . 

1 2 




S YIX = 


£ 


/log w c \ 2 


/ 


= log w + fAw c 1 


(43) 


i=l 


\ ~h 


/ 13 


-1 



where the final term in Eq. (43) represents the mean scatter about the 
A 

graph of Y = AX. 

For example, at 0000 GMT, 21 March, 1971, the standard error at Way- 

cross, Ga . , was 0.1131. Therefore, .1131 = log w + |A W C I , hence it 

w 

follows that w 4 * | A w c I = 1.298. This means that the mean residual ratio, 

w 

after regression, for this case is | A w c * = 30%. Analysis of the actual 

residuals for the soundings shown indicates that the majority of this error- 
ratio occurs at the highest altitudes recorded in the sounding where w is 
the smallest and where actual instrument accuracy is doubtful. 
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TABLE III 



Values of A. and related statistics 
for the data-sample for Waycross, Ga . 



DATE 
DAY -GMT 


X -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

X -VALUE 


CORRELATION 

COEFFICIENT 

Rl 2 Rl 2 


STANDARD 

DEVIATION 

cry 


STANDARD 

ERROR 

Syix 


16-0000 


2.9414 


0.2585 


.953 


.908 


0.5861 


0.1837 


16 - 1200 * 


1.5298 


0.4277 


.704 


.496 


0.4130 


0.3043 


17-0000 


4.5798 


0.1950 


.988 


.976 


0.8809 


0.1387 


17-1200 


2.7855 


0.1246 


.987 


.974 


0.5373 


0.0888 


18-0000 


2.2575 


0.1187 


.983 


.966 


0.4377 


0.0846 


18-1200 


2.6057 


0.1371 


.983 


. 966 


0.5057 


0.0978 


19-0000 


2.3259 


0.3495 


.879 


.773 


0.5038 


0.2491 


19-1200 


1.8176 


0.2879 


.868 


.755 


0.3981 


0.2049 


20-0000 


3.6279 


0.2036 


.980 


.960 


0.7030 


0.1447 


20-1200 


2.5391 


0.2308 


.950 


.903 


0.5083 


0.1643 


21-0000 


3.0887 


0.1589 


.983 


.966 


0.5976 


0.1131 


21-1200 


2.4250 


0.1087 


.987 


.974 


0.4680 


0.0775 


22 - 0000 * 


0.9815 


0.3252 


.642 


.412 


0.2911 


0.2316 


22-1200 


3.2945 


0.2850 


.955 


.912 


0.6572 


0.2031 


23-0000 


4.3912 


0.2107 


.985 


.970 


0.8473 


0.1499 


23-1200 


4.5008 


0.2425 


.982 


. 964 


0.8713 


0.1724 


TIME MEAN 


2.8557 






.869 






STD. DEV. 
X ) 


1.0401 













^Corresponds to < -85. 
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TABLE IV 



Values of A and related statistics 
for the data-sample for Charleston, S. C. 



DATE 
DAY -GMT 


X -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

X -VALUE 


CORRELATION 

COEFFICIENT 

Rl2 Rl2 


STANDARD 

deviation 

a Y 


STANDARD 

ERROR 

SY|X 


16-0000 


2.3199 


0.1174 


.984 


.968 


0.4485 


0.0863 


16-1200 


4.33721 


0.2933 


.972 


.945 


0.8494 


0.2088 


17-0000 


3.9920 


0.1322 


.993 


.986 


0.7646 


0.0941 


17-1200 


3.8303 


0.1629 


.989 


.978 


0.7379 


0.1160 


18-0000 


2.4849 


0.1724 


.970 


.941 


0.4883 


0.1229 


18-1200 


2.7495 


0.1236 


.987 


.974 


0.5319 


0.0883 


19-0000 


2.9405 


0.1308 


.988 


.976 


0.5684 


0.0934 


19-1200* 


1.5187 


0.4050 


.721 


.520 


0.4013 


0.2886 


20-0000 


4.5054 


0.2364 


.983 


.966 


0.8711 


0.1680 


20-1200 


2.7135 


0.2153 


.961 


.924 


0.5369 


0.1532 


21-0000 


4.7245 


0.1834 


.990 


.980 


0.9077 


0.1305 


21-1200 


2.4194 


0.1502 


.976 


.953 


0.4727 


0.1071 


22-0000 


MISSING 


DATA 










22-1200 


2.5801 


0.2807 


.931 


.867 


0.5281 


0.2001 


23-0000 


4.4411 


0.2290 


.983 


.966 


0.8597 


0.1631 


23-1200 


3.5360 


0.1560 


.988 


.976 


0.6809 


0.1110 


TIME MEAN 


3.2729 






.928 






STD. DEV. 

V A > 


0.9832 













^Corresponds to R-^ < .85. 
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TABLE V 



Values of A and related statistics 
for the data-sample for Lake Charles, La. 



DATE A -VALUE 

DAY -GMT CALCULATED 


STANDARD 
ERROR OF 

A -VALUE 


CORRELATION 

COEFFICIENT 

R 12 R 12 


STANDARD 

DEVIATION 

<7 Y 


STANDARD 

ERROR 

S Y 1 X 


16-0000 


2.9224 


0.1461 


.984 


.968 


0.5653 


0.1041 


16-1200 


4.0954 


0.1065 


.996 


.992 


0.7837 


0.0759 


17-0000 


4.5898 


0.1810 


.990 


.980 


0.8835 


0.1291 


17-1200 


3.5190 


0.2167 


.976 


.953 


0.6878 


0.1547 


18-0000 


4.0538 


0.1916 


.986 


.972 


0.7843 


0.1367 


18-1200 


3.8884 


0.1964 


.984 


.968 


0.7528 


0.1399 


19-0000 


2.7501 


0.0960 


.992 


.984 


0.5273 


0.0684 


19-1200 


3.0506 


0.2307 


.965 


.931 


0.6030 


0.1646 


20-0000 


2.4835 


0.0959 


.990 


.980 


0.4782 


0.0684 


20-1200 


2.8351 


0.1944 


.971 


.943 


0.5575 


0.1388 


21-0000* 


1.1751 


0.2333 


.813 


.661 


0.2757 


0.1665 


21-1200 


1.6059 


0.2349 


.885 


.783 


0.3462 


0.1676 


22-0000 


4.3386 


0.1766 


.992 


.984 


0.7777 


0.1041 


22-1200 


5.2771 


0.1789 


.994 


.988 


0.8705 


0.1016 


23-0000 


4.4083 


0.2264 


.983 


.966 


0.8532 


0.1612 


23-1200 


4.5449 


0.3695 


.963 


.927 


0.8417 


0.2375 


TIME MEAN 
STD. DEV. 

M A ) 


3.4711 

1.1352 






.936 






■^Corresponds 


to ^ 


.85. 
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TABLE VI 



Values of A and related statistics 





for 


the data-sample for 


El Paso, 


Texas . 




DATE 
DAY -GMT 


A -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

A -VALUE 


CORRELATION 

COEFFICIENT 

R 12 R 12 


STANDARD 

DEVIATION 

a Y 


STANDARD 

ERROR 

S YIX 


16-0000 


3.5147 


0.3078 


.964 


.929 


0.5391 


0.1509 


16-1200 


3.3290 


0.2334 


.976 


.953 


0.5050 


0.1147 


17-0000 


3.4696 


0.1695 


.988 


.976 


0.5195 


0.0832 


17-1200 


3.2455 


0.3847 


.936 


.876 


0.5123 


0.1886 


18-0000 


3.7536 


0.2578 


.997 


.994 


0.5658 


0.1260 


18-1200 


3.3012 


0.1721 


.987 


.970 


0.4938 


0.0843 


19-0000 


3.3597 


0.1046 


.995 


.990 


0.4990 


0.0513 


19-1200 


2.4246 


0.2875 


.936 


.876 


0.3841 


0.1414 


20-0000 


2.5997 


0.2698 


.950 


.902 


0.4057 


0.1327 


20-1200* 


0.6696 


0.1802 


.778 


.605 


0.1311 


0.0868 


21-0000 


2.6959 


0.3116 


.939 


.882 


0.4232 


0.1524 


21-1200 


3.1878 


0.1860 


.983 


.966 


0.4787 


0.0911 


22-0000 


3.3213 


0.2029 


.982 


.964 


0.4995 


0.0994 


22-1200 


2.4686 


0.3511 


.912 


.832 


0.4003 


0.1722 


23-0000 


2.3344 


0.1755 


.973 


.947 


0.3545 


0.0860 


23-1200* 


0.7583 


0.4089 


.526 


.277 


0.2192 


0.1965 


TIME MEAN 

STD. DEV. 
ff t ( A ) 


2.7771 

0.9147 






.871 







Corresponds to < *85. 



TABLE VII 



Values of A and related statistics 
for the data-sample for San Diego, Calif. 



DATE 
DAY -GMT 


A -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

A -VALUE 


correlation 

COEFFICIENT 

2 

R 12 R 12 


STANDARD 

DEVIATION 

a Y 


STANDARD 

ERROR 

S YI X 


16-0000 


3.5187 


0.0897 


. 996 


.992 


0.6713 


0.0637 


16-1200 


3.8860 


0.1442 


.991 


.982 


0.7444 


0.1024 


17-0000 


3.4676 


0.0874 


.996 


.992 


0.6606 


0.0620 


17-1200 


3.2899 


0.1467 


.988 


.976 


0.5932 


0.0943 


18-0000 


3.9886 


0.1914 


.987 


.974 


0.7205 


0.1230 


18-1200 


4.0171 


0.1426 


.992 


.984 


0.7684 


0.1012 


19-0000 


3.6748 


0.1190 


.993 


.986 


0.7023 


0.0845 


19-1200 


2.4739 


0.1570 


.975 


.951 


0.4817 


0.1115 


20-0000* 


1.9672 


0.3994 


.807 


.651 


0.4627 


0.2836 


20-1200 


3.9124 


0.1845 


.986 


.972 


0.7529 


0.1309 


21-0000 


3.3273 


0.2063 


.976 


.953 


0.6470 


0.1465 


21-1200 


3.8186 


0.1946 


.984 


.968 


0.7370 


0.1382 


22-0000 


2.6873 


0.2288 


.956 


.914 


0.5338 


0.1626 


22-1200 


2.3642 


0.2423 


.938 


.880 


0.4786 


0.1721 


23-0000 


2.1203 


0.2521 


.919 


.845 


0.4378 


0.1790 


23-1200 


2.8397 


0.1020 


.992 


.984 


0.5434 


0.0725 


TIME MEAN 


3.2096 






.938 







STD. DEV. 0.7012 
a t ( A. ) 

Corresponds to R^2 < -85. 
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TABLE VIII 



Values of A and related statistics 
for the data-sample for Caribou, Maine. 



DATE 
DAY -GMT 


A -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

A -VALUE 


CORRELATION 

COEFFICIENT 

R 12 R 12 


standard 

DEVIATION 

°X 


STANDARD 

ERROR 

S Y IX 


16-0000* 


1.3979 


0.2714 


.841 


.707 


0.2937 


0.1661 


16-1200 


4.3499 


0.0839 


.998 


.996 


0.7705 


0.0534 


17-0000 


3.5403 


0.1644 


.988 


.976 


0.6588 


0. 1047 


17-1200 


2.9646 


0.1566 


.984 


.968 


0.5350 


0.1003 


18-0000 


3.1736 


0.2673 


.960 


.922 


0.5869 


0.1711 


18-1200* 


1.5775 


0.4351 


.723 


.523 


0.3878 


0.2788 


19-0000 


4.6004 


0.1396 


.995 


.990 


0.8453 


0.0884 


19-1200 


3.7880 


0.2120 


.982 


.964 


0.6859 


0.1359 


20-0000* 


1.0878 


0.3522 


. 666 


.444 


0.2901 


0.2254 


20-1200 


2.7297 


0.2756 


.957 


.916 


0.4660 


0.1424 


21-0000 


5.0908 


0.3520 


.979 


.958 


0.9154 


0.1960 


21-1200 


5.1735 


0.2328 


.989 


.978 


0.9479 


0.1461 


22-0000 


3.7348 


0.0533 


.999 


.998 


0.6789 


0.0335 


22-1200 


4.5463 


0.1532 


.994 


.988 


0.8320 


0.0965 


23-0000 


6.5325 


0.5259 


.975 


.951 


1.1656 


0.2745 


23-1200 


4.5879 


0.1398 


.996 


.992 


0.8014 


0.0769 


TIME MEAN 

STD. DEV. 
ff t( A ) 


3.6797 

1.4896 






.892 







Corresponds to R-^ < .85. 
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TABLE IX 



Values of X and related statistics 
for the data-sample for Maniwaki, Que. 



DATE X -VALUE 

DAY -GMT CALCULATED 


STANDARD 
ERROR OF 

X -VALUE 


CORRELATION 

COEFFICIENT 

r 12 r 12 


STANDARD 

DEVIATION 

°Y 


STANDARD 

ERROR 

S YIX 


16-0000 


5.5686 


0.32 28 


.981 


.962 


1.0019 


0.2053 


16-1200 


3.0656 


0.2679 


.971 


.943 


0.4894 


0.1246 


17-0000 


3.6455 


0.1713 


.987 


.974 


0.6561 


0.1097 


17-1200 


3.2185 


0.2048 


.977 


.955 


0.5867 


0.1315 


18-0000 


3.2674 


0.2213 


.976 


.953 


0.6138 


0.1405 


18-1200 


2.8017 


0.1775 


.977 


.955 


0.5597 


0.1249 


19-0000 


4.2382 


0.1084 


.996 


.992 


0.8075 


0.0770 


19-1200* 


1.1753 


0.3235 


.724 


.524 


0.2892 


0.2077 


20-0000 


1.8514 


0.4214 


.874 


.764 


0.3945 


0.2075 


20-1200 


3.0911 


0.1692 


.987 


.974 


0.5533 


0.0945 


21-0000 


3.2247 


0.0641 


.999 


.998 


0.5703 


0.0320 


21-1200 


2.8132 


0.1215 


.991 


.982 


0.5029 


0.0714 


22-0000 


2.3573 


0.2628 


.933 


.870 


0.4481 


0.1680 


22-1200 


2 . 6480 


0.1180 


.989 


.978 


0.4888 


0.0747 


23-0000 


3.4791 


0.1644 


.989 


.978 


0.6695 


0.1038 


23-1200 


3.8057 


0.3301 


.971 


.943 


0.6607 


0.1670 


TIME MEAN 


3.1407 






.922 






STD. DEV. 
o t ( A ) 


0.9849 












^Corresponds 


to R^2 < 


.85. 
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TABLE X 



Values of X and related statistics 
for the data-sample for Sault-S te-Marie , Mich. 



DATE 
DAY -GMT 


A -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

A -VALUE 


CORRELATION 

COEFFICIENT 

R 12 r 12 


STANDARD 

DEVIATION 

° Y 


STANDARD 

ERROR 

S YIX 


16-0000 


2.4004 


0.3536 


.932 


.869 


0.4486 


0.1741 


16-1200 


2.4229 


0.1205 


.986 


.972 


0.4355 


0.0770 


17-0000 


4.3460 


0.1218 


.995 


.990 


0.7756 


0.0780 


17-1200 


3.6100 


0.1099 


.995 


.990 


0.6639 


0.0697 


18-0000 


3.7061 


0.1273 


.994 


.988 


0.6828 


0.0807 


18-1200 


3.5366 


0.1136 


.994 


.988 


0.6515 


0.0721 


19-0000 


2.0358 


0.2896 


.904 


.817 


0.3889 


0.1733 


19-1200 


2.2482 


0.2289 


.943 


.889 


0.4225 


0.1463 


20-0000 


2.8100 


0.1758 


.979 


.958 


0.5218 


0.1108 


20-1200 


3.2512 


0.0385 


.999 


.998 


0.5917 


0.0242 


21-0000 


3.3835 


0.1394 


.991 


.982 


0.6221 


0.0879 


21-1200 


2.1833 


0.2288 


.945 


.893 


0.4217 


0. 1446 


22-0000 


3.0157 


0.1823 


.981 


.962 


0.5614 


0.1153 


22-1200 


3.3318 


0.1228 


.993 


.986 


0.6384 


0.0775 


23-0000 


3.8431 


0.1015 


.997 


.994 


0.7333 


0.0640 


23-1200 


3.7412 


0.1077 


.996 


.992 


0.7496 


0.0680 


TIME MEAN 


3.1166 






.954 







STD. DEV. 0.6940 

M A > 
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TABLE XI 



Values of A and related statistics 
for the data-sample for Bismarck, N. D. 



DATE 
DAY -GMT 


STANDARD 
ERROR OF 

A -VALUE 

CALCULATED A. -VALUE 


CORRELATION 

COEFFICIENT 

R 12 R 12 


STANDARD 

DEVIATION 

a Y 


STANDARD 

ERROR 

S Yl X 


16-0000 


2.6549 


0.2499 


.955 


.912 


0.5027 


0.1564 


16-1200 


3.4759 


0.2024 


.984 


.968 


0.6213 


0.1180 


17-0000 


6.5097 


0.1724 


.996 


.992 


1.1819 


0.1080 


17-1200 


NO DATA 


AVAILABLE 










18-0000 


NO DATA 


AVAILABLE 










18-1200 


3.6856 


0.2127 


.984 


.968 


0.6601 


0.1243 


19-0000 


4.4516 


0.2273 


.986 


.972 


0.8179 


0.1427 


19-1200 


4.5117 


0.1570 


.995 


.990 


0.8724 


0.0955 


20-0000 


3.4355 


0.1859 


.983 


.966 


0.6154 


0.1180 


20-1200 


2.0076 


0.2615 


.918 


.843 


0.3888 


0.1611 


21-0000 


4.9089 


0.1923 


.992 


.984 


0.8086 


0.1088 


21-1200 


4.5735 


0.1341 


.996 


.992 


0.7678 


0.0743 


22-0000 


5.0893 


0.2730 


.987 


.974 


0.9794 


0.1640 


22-1200 


3 . 8804 


0.0609 


.999 


.998 


0.7341 


0.0382 


23-0000 


3.6326 


0.3285 


.954 


.910 


0.6710 


0.2088 


23-1200 


2.1195 


0.1005 


.987 


.974 


0.3789 


0.0639 


TIME MEAN 


3.9240 






.960 







STD. DEV.’ 1.2147 
* t ( * ) 
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TABLE XII 



Values of A and related statistics 
for the data-sample for Glasgow, Mont. 



DATE 
DAY -GMT 


A -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

X -VALUE 


CORRELATION 

COEFFICIENT 

R 12 R 12 


STANDARD 

DEVIATION 

°Y 


STANDARD 

ERROR 

S YI X 


16-0000 


3.7842 


0.1697 


.990 


.980 


0.6406 


0.0943 


16-1200 


6.6178 


0.3483 


.986 


.972 


1.1253 


0.1938 


17-0000 


5.6933 


0.1384 


.997 


.994 


0.9578 


0.0770 


17-1200 


6.0083 


0.2198 


.993 


.986 


1.0154 


0.1224 


18-0000 


6.0189 


0.1978 


.995 


.990 


1.0630 


0.1100 


18-1200 


5.1013 


0.1940 


.993 


.986 


0.8650 


0.1083 


19-0000 


6.0881 


0.1843 


.995 


.990 


1.0284 


0.1028 


19-1200 


2.6928 


0.2074 


.972 


.945 


0 . 4666 


0.1158. 


20-0000 


3.3309 


0.2818 


.963 


.927 


0.5644 


0.1593 


20-1200 


2.2316 


0.1834 


.965 


.931 


0.3754 


0.1031 


21-0000 


5.3035 


0.2107 


.991 


.982 


0.8700 


0.1187 


21-1200 


3.1120 


0.3134 


.953 


.908 


0.5487 


0.1746 


22-0000 


4.4330 


0.1457 


.995 


.990 


0.7500 


0.0813 


22-1200 


3.3728 


0.1415 


.991 


.982 


0.5573 


0.0802 


23-0000* 


0.6804 


0.3303 


.614 


.377 


0.1954 


0.1649 


23-1200* 


0.8355 


0.2358 


.746 


.557 


0.1885 


0.1316 


TIME MEAN 

STD. DEV. 
°t( X ) 


4.0815 

1.8749 






.906 







Corresponds to R-^ < *85. 
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TABLE XIII 



Values of X and related statistics 
for the data-sample for Spokane, Wash. 



DATE 
DAY -GMT 


X -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

X -VALUE 


CORRELATION 

COEFFICIENT 

2 

R 12 R 12 


STANDARD 

DEVIATION 

°Y 


STANDARD 

ERROR 

S YIX 


16-0000 


6.1242 


0.2171 


.994 


.988 


1.0350 


0.1209 


16-1200 


5.3894 


0.0292 


.999 


.998 


0.9060 


0.0163 


17-0000 


6.2015 


0.2579 


.992 


.984 


1.0971 


0.1432 


17-1200 


5.4573 


0.0661 


.999 


.998 


0.9586 


0.0367 


18-0000 


4.1476 


0.2953 


.973 


.947 


0.6969 


0.1673 


18-1200 


3.1086 


0.2754 


.959 


.920 


0.5305 


0.1562 


19-0000 


3.3134 


0.1454 


.990 


.980 


0.5484 


0.0825 


19-1200 


2.5267 


0.1350 


.985 


.970 


0.4199 


0.0766 


20-0000 


2.6138 


0.1618 


.980 


.960 


0.4352 


0.0914 


20-1200 


4.5322 


0.4601 


.948 


.899 


0.7780 


0.2593 


21-0000 


4.2112 


0.3912 


.956 


.914 


0.7166 


0.2204 


21-1200 


3.6604 


0.2030 


.984 


.968 


0.6057 


0.1145 


22-0000* 


1.6059 


0.3356 


.822 


.676 


0.3184 


0.1894 


22-1200* 


1.0266 


. 0.2085 


.829 


.687 


0.2015 


0.1176 


23-0000 


5.3603 


0.2307 


.990 


.980 


0.8799 


0.1299 


23-1200 


4.0725 


0.4338 


.968 


.937 


0.5157 


0.1406 ' 


TIME MEAN 


3.9594 






.925 







STD. DEV. 1.5414 
ff t( * > 

•'Corresponds to R-^ < -85. 
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TABLE XIV 



Values of A and related statistics 
for the data-sample for Quillayute, Wash. 



DATE 
DAY -GMT 


A -VALUE 
CALCULATED 


STANDARD 
ERROR OF 

A -VALUE 


CORRELATION 

COEFFICIENT 

r 12 R 12 


STANDARD 

DEVIATION 

°Y 


STANDARD 

ERROR 

S YIX 


16-0000 


5.5314 


0.1401 


.996 


.992 


1.0899 


0.0992 


16-1200 


5.5505 


0.1507 


.996 


.992 


1.1368 


0.1065 


17-0000 


4.3589 


0.2506 


.979 


.958 


0.8494 


0.1789 


17-1200* 


1.4652 


0.3352 


.771 


.594 


0.3626 


0.2394 


18-0000* 


1.7747 


0.3777 


.793 


.629 


0.4275 


0.2701 


18-1200 


1.6247 


0.2586 


.867 


.752 


0.3580 


0.1850 


19-0000 


2.0580 


0.1481 


.968 


.937 


0.4061 


0.1058 


19-1200 


2.3200 


0.1796 


.963 


.927 


0.4596 


0.1282 


20-0000 


2.5319 


0.1793 


.969 


.939 


0.4977 


0.1278 


20-1200 


3.6523 


0.3188 


.954 


.910 


0.7291 


0.2272 


21-0000* 


1.6818 


0.2995 


.842 


.709 


0.3805 


0.2133 


21-1200 


1.4532 


0.3086 


.794 


.630 


0.3485 


0.2199 


22-0000 


1.7184 


0.2594 


.878 


.771 


0.3722 


0.1846 


22-1200 


2.0152 


0.1758 


.982 


.964 


0.4858 


0.1019 


23-0000 


4.1114 


0.1946 


.986 


.972 


0.7917 


0.1382 


23-1200 


3.9581 


0.3402 


.972 


.945 


0.7002 


0.1755 


TIME MEAN 


2.8629 






.851 







STD. DEV. 1.5414 
ff t( A ) 

Corresponds to R 12 < .85. 
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In table III, a low correlation of 0.642 occurs at 22 March, 0000GMT, 
which results in a relatively high standard error Sy|x- Note in Table III 
that most of the profiles have with associated small standard 

errors. There were two cases in Table III with low correlation coefficients, 
specifically 0.704 and 0.642. This pattern of correlation coefficients 
Rl2 ^ 0.85 except for a possible set of one to three profiles for which 
R ^ fell below this value, continues throughout all of the 12 stations 
considered here as shown in Tables III through XIV. Specific instances in 
these tables were the R-j^-value value drops below 0.85 have been marked 
by an asterisk alongside of the date. In all of these cases the parameter 
X is likewise smaller than normal in view of Eq. (20) 



All of these special cases where R 12 510.85 are characterized by a 
standard deviation a y lower than typical for the station. Such cases of 
small a y-values represented cases where statistically there was insufficient 
vertical variation in the moisture range to ascribe it successfully to the 
power- prof ile parameterization. 

An example of a sma Her- than-norma 1 correlation coefficient R-^ and 
its correspondingly small X has been extracted from the regression print 
out in the form of Table XV. This table lists the statistical results 
from Waycross, Ga . , at 0000GMT, 22 March, 1971. This power-profile 
analysis resulted in R-^ = 0.642 and X = 0.9815, which was considerably 
lower than the mean X -value. In Table XV, the observed moisture profile 
is indicated in column 2 as Y of Eq. (7), and in column 3 is shown the 



corresponding Y of Eq. (40). The residual shown in column 4 results from 
A A 

forming Y-Y from Eq. (41). As expected, Y began as a maximum positive 
value at the surface and decreased to 0.000 at 500 mb and becomes negative 




(44) 



A 
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TABLE XV 



An example for Waycross, Ga . (0000 GMT, 22 Mar., 1971) 
showing a condition leading to an ill-defined A -profile. 



PRESSURE (mb) 


Y (obs) 


A 

Y 


RESIDUAL 


1014 


0.4660 


0.3014 


0.1646 


1000 


0.3179 


0.2955 


0.0224 


950 


0.2719 


0.2736 


7 O . 0017 


900 


0.2478 


0.2505 


-0.0028 


850 


0.2020 


0.2262 


-0.0241 


800 


0.1024 


0.2003 


-0.0979 


750 


-0.0060 


0.1728 


-0.1789 


700 


-0.1507 


0.1434 


-0.2942 


650 


-0.1579 


0.1118 


-0.2716 


600 


0.1296 


0.0777 


0.0519 


550 


-0.0082 


0.0406 


-0.0488 


500 


0.0000 


0.0000 


0.0000 


450 


-0.3650 


-0.0449 


-0.3201 


400 


-0.6949 


-0.0951 


-0.5998 
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1000 

950 

900 

850 

800 

750 

700 

650 

600 

550 

500 

450 

400 



TABLE XVI 



Environmental parameters used in the 
calculation of Table XV. 



TEMP (C) 


RELATIVE 
HUMIDITY, % 


MIXING 
RATIO (g/kg) 


16.8 


33 


3.8801 


15.6 


25 


2.7590 


12.0 


27 


2.4816 


8.3 


31 


2.3477 


6.9 


29 


2.1130 


4.3 


26 


1.6800 


1.6 


23 


1.3087 


- .5 


18 


.9378 


- 2.9 


20 


.9186 


- 6.3 


48 


1.7885 


-10.9 


48 


1.3021 


-16.9 


77 


1.3270 


-19.7 


39 


.5726 


-27.2 


34 


.2678 
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above 500 mb. The observed profile, Y (p) , began also as a positive value 
at the surface, but becomes negative in the layer 750 to 600 mbs indicating 
that the actual mixing ratio profile had values less than W 500 in that layer. 

Table XVI lists the actual a tmospheric- sounding parameters which resulted 
in the values shown in Table XV. Note that the reported relative humidities 
are low in the 750 to 650 mb layer and high in the layer 600 to 450 mb. 

Since by Eq. (9), w is a function of relative humidity as well as of 
temperature and pressure, relatively low w-va lues exis t from 750 to 650 mb 
and high w-values in the 600- to 450 mb- range. This, at least locally, is 
inconsistent with the assumption on which the lower-fit is based and, as 
calculated, gives a low value of R-j^. Note that the power-fit need not 
have been inconsistent if the temperature lapse was large enough through 
the 600 - 450 mb layer 5 which in this instance was not the case. 

As just noted, the second major factor which controls the consistency 
of the profile-fit at level p is the temperature T(p). It can also be 
shown, by example, that a temperature inversion or an isothermal layer 
may be (depending upon relative humidity) associated with w values which 
do not decrease in height, resulting in the low correlation coefficients 
Rf 2 . In Tables III to XIV, all cases where R-j^ - 0.85 which thus led to 
an ill-defined (i.e., too small ) value of X are denoted by an asterisk. 

The explanation of all of these special cases is similar to that advanced 
for the Waycross, Ga . case in Tables XV, XVI. 



Note finally that it is convenient to list the time-mean X and R-j^ 
for each station listed in Tables III,..., XIV. In addition it is desir- 
able to compute the temporal standard deviation of the X 1 s in each 
table according to 



2 



a t ( A ) = ( A 1 - A ) 2 +. . . + ( 




( 45 ) 



15 
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B. REDUCTION OF VARIANCE OF A 



Having obtained the sixteen A -values per station the next portion 
of the investigation involved determining the feasibility of reducing the 
variance of A for each station over the period of time considered using 
other sounding parameters. The synoptic variables introduced into the 
computation as independent variables permitted utilization of the stepwise 
multiple regression aspect of the BIMED 02R program where A was the 
dependent variable. In addition to the 16 values of A per station, 16 
values of each of the following were read into the program as indepen- 
dent variables: 

surface temperature defined as 
surface mixing ratio defined as X 2 
850 mb temperature defined as X^ 

850 mb mixing ratio defined as X^ 

700 mb temperature defined as X^ 

500 mb temperature defined as X^. 

The above variables were used in the multiple regression as inputs with 
the computed A as the dependent variable in each case. The desired 
solution equation would then be of the form of Eq. (32) where « A , 
and the X^ are the respective independent variables. The stepwise 
regression program BIMED 02R was then utilized to derive the best-fit 
form of 

A = A Q + AjX x + A 2 X 2 +. . .+A 6 X 6 (46) 

Among the factors which might be considered as useful in reducing the 
temporal variance of A are variations in the synoptic patterns, and those 
which may reflect purely local conditions. Among the latter effects 
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(e.g., Tables III,..., XIV), station location and elevation should be 

considered, as well as proximity to large water bodies and lakes. In 

the sounding variables X , , . ,X , the lower-level parameters, were con- 
i' 6 

sidered to be primarily representative of local effects, whereas the 
higher-level variables were presumed to be indicative of the larger-scale 
synoptic time-scales, although there is undoubtedly some interaction 
between the low- and high-level variables. 

The temporal standard deviations cr^_( A ) listed at the base of 
Tables III,..., XIV resulting from Eq. (45) ranged from minimal values 
of ^ t ( A ) = 0.694 at Sault Ste Marie, and o A ) = 0.701 at 
San Diego to a maximum o t ( A ) = 1.875 at Glasgow, Montana. In general 
it can be seen that stations close to large bodies of water (Charleston, 
San Diego, Maniwaki, and Sault Ste Marie) have relatively small standard 
deviations whereas stations away from bodies of water have larger standard 
deviations (Bismarck, Glasgow, and Spokane). The cases of large values 
of <7 t ( A ) seem to afford considerable opportunity for utilization 
of the regression procedure for the reduction of the A -variance. How- 
ever, it was decided to test the diagnostic usefulness of Eq. (46) for all 
twelve stations, regardless of the o ^_( A ) value. 

A preliminary test run was conducted using the BIMED 02R program with 
the variables defined previously in the section. During this run all 
independent variables were allowed to enter freely based solely on their 
partial correlation values. The program was allowed to run until the 
requirement of 2% added fractional explained variance after introduction 
of the kth variable could not be met. 

Two facts became evident from this preliminary test: (1) the indepen- 

dent variable with the highest initial correlation coefficient in almost 
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all cases was X^, and (2) the maximum number of significant steps necessary 
to ensure the required minimum added explained variance at each step was 
four. At steps five and six of the preliminary run the added explained 
variance upon predictor-entry was too small. Since a low-level moisture 
predictor was desired to be compatible with the variations of the X - 
profile parameter (cf., the discussion for Waycross at 0000 GMT, 22 March, 
1971), it was decided to utilize the control-delete option of BIMED 02R 
to force in X^ as the first predictor at each station. The second option 
used was that of limiting the number of predictors to four in each 
regression analysis. The results of the sets of regressions over the 12 
stations are shown in Table XVII. 

Equation (42) gives the relationship between the standard error of 
X , after the four-step multiple correlation has been determined 

S A | X = (L " V A > (47) 

From Table XVII it can be seen that when a^( A ) is relatively small, 

the explained fraction of this statistic is relatively small. The result, 

then, of application of the multiple regression technique is to reduce 

the standard deviation O ( A ) of all stations, but to reduce the 

stations with the largest ^ ( A ) the most. 

San Diego is a good example of a station having a small A )> 

a t ( A ) = 0.7012, which means, according to the above result, a 
relatively small reduction, specifically to S>^|^=0.5686. Correspondingly, 
this small additional reduction is reflected in the relatively low multiple 
correlation coefficient at San Diego of R m = .5852. At the other extreme 
is Glasgow having a temporal standard deviation o A ) = 1.8749. By 
the multiple regression procedure a standard error of estimate of 0.8745 
results. Again this is consistent with the high multiple correlation 
coefficient, R m = 0.8846. 
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Stepwise regression parameters defining A in terms 
sounding parameters for the data period (16 soundings). 
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1394 




Figure 2: 



Histogram depicting number of stations within certain 
standard error limits. 
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The result of the multiple regression for each station is to give a re- 
duced standard error of A . A histogram depicting the 12 standard error 
values is shown in Figure 2, on a scale of S from 0.50 to 0.95. Note 

that most of the standard errors after regression fall into the range 
0.5 to 0.65. 

The righthand column in Table XVII lists the order of entry of the 
independent variables into the stepwise regression and is based on the 
value of the partial correlation of the respective variables as explained 
in Section III (B) . From Table XVII it may be noted that in seven of the 
12 stations the second independent variable added to the regression was X^. 
Also in seven of the 12 stations the third independent variable added was 
X 2 . This led to the idea of possible classification of stations according 
to a forced order of entry, an option of BIMED 02R, and will be discussed 
in the next sub-section. 

C. STATION CLASSIFICATION ACCORDING TO THE PREDICTOR ENTRY 

Table XVIII lists an attempt to classify stations by the relative 
priority of the independent variables to the program specified in the order 
X 4 , X^, X^, and the remaining predictor in a 4-step selection left as free 
choice according to the selection criterion. As can be seen by comparing 
Tables XVII and XVIII there was no change in the multiple correlation 
coefficient at 10 of the 12 stations as a result of specifying the 4 - 6-2 
input order. This is because in almost all cases the variables X 4 , Xg, and 
X 2 were in the original run (Table XVII) even though not necessarily in this 
specific order. Table XVIII lists the "single free choice" variable under 
order of input as the last independent variable in the list for each station. 
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Stepwise regression parameters defining A in terms of sounding parameter 
in a specified order for the data period (16 soundings). 
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MEAN ff t ( A ) = 1.1394 



It may be noted that there are four stations with surface temperature as 
the last independent variable (Charleston, Maniwaki, Bismarck, and 
Spokane), four with 500 mb- tempera ture (Waycross, Lake Charles El Paso, 
and Caribou), and four with 850 mb- temperature (San Diego, Sault Ste 
Marie, Glasgow, And Quillayute) . There might, then, be some way of 
classifying or typing stations on the basis of the free-choice predictor 
selected at step 4. However, with the small station set available it 
was not clear that this could be done. 

Finally, it is to be emphasized that Table XVII indicates that the 
observed A -values tend to adjust to (a) the synoptic-scale time 
variations in the sounding and to (b) the local-effect variations. The 
latter effect is harder to specify so that one achieves some limiting 
accuracy in the expected value resulting from the regression. The magnitude 
of this accuracy is suggested by the standard error which tends to show a 
cluster of values in the range 0.50 to 0.65 (Fig. 2). 
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V. PRECIPITABLE WATER VAPOR TESTS 



As noted in the introduction, two verification techniques were to be 
used to evaluate the results of this investigation. The first of these, 
namely verification of the moisture profile, has been discussed in detail 
in Section IV from a statistical viewpoint. In this section, an analysis 
of the verification of the moisture-parameterization technique from the 
viewpoint of precipitable water vapor is presented. 

The formulas used to calculate the precipitable water are shown as 
Eqs. (12) and (13). Equation 12 was utilized to obtain precipitable water 
based upon the moisture profile of the data for the given data sample. It 
is the precipitable water from the observed sounding, and hence will be 
referred to as observed precipitable water, W Q . Equation 13 is used to 
calculate the precipitable water utilizing the X -values already obtained. 
As can be seen from Eq. (13), it is based on only three variables: 500-mb 

mixing-ratio, surface pressure, and X . There are, however, two possible 
methods since X -values were obtained by the profile-parameter method 
and, secondly, by the multiple regression technique. Precipitable water 
calculated using the profile parameter X -values, one for each sounding, 
will be designated as W cs and those relatively smooth values obtained by 
using the multiple regression X -values denoted by W . 

For this phase of the investigation, four stations were chosen primarily 
on the basis of their varying geographical locations: Waycross, Ga . , in 

the southeastern United States; San Diego in the southwest; Spokane, Wash., 
in the northwest; and Sault Ste Marie, Mich., in the northeast. Results 
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of the precipitable water calculations are listed in Tables XIX through 
XXII and are plotted in Figures 3 through 6. 

As can be seen from Tables XIX, ..., XXII or from Figures 3,..., 6, the 
W cs -value is in almost all cases within approximately 0.05 cms of the W Q 
value. Needless to say, with such good verification, the trend of W cs 
is the same as the observed trend. It is obvious from the results, then, 
that the use of the power-profile technique can be expected to yield 
precipitable water values which are very representative of the W Q values, 
i.e. within 0.05 cms (cf., Tables XIX, . . . ,XXII) . 

The values of W cm are also presented in Tables XIX, . . . ,XXII and in 
Figures 3,..., 6. As can be seen errors do exist in the W cm values, 
assuming W Q is the correct value. These errors can at times become large, 
i.e. as much as 1007 o . In general, however, the W cm values which result 
from using X -values obtained by multiple regression follow the time 
trend of W . Large errors in W cm are then the exception rather than the 
rule, and as can be seen in Figures 3 ,..., 6 , occur at only one or two 
sounding times per station. The soundings which give the largest W cm 
errors are, in fact, just those soundings which had the lower values of 
X and of the correlation coefficients R^2 i n connection with the W cs 
computations . 

The 0000 GMT, 22 March, 1971 sounding at Waycross, Ga . will be used 
to illustrate the point of a large W C m error. As can be seen in Table III, 
this sounding gave a low correlation coefficient value, 0.642, and a X - 
value of 0.9815 where the time-mean over the 16 soundings for that station 
was 2.8557. Figure 3 shows that the use of the X -value from the power- 
fit profile, X = 0.9815, in Eq. (13) gives an accurate verification of 
W Q . On the other hand the value of X resulting from the multiple 
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TABLE XIX 



Table of precipitable water vapor values (W q> W cs , W cm ) 

for Waycross, Ga . 



DAY 

DAY - GMT 


OBSERVED 

W o 


PROFILE 

W 

w cs 


MULT. REG. 
W cm 


16 - 0000 


3.7190 


3.8704 


3.0169 


16 - 1200 


1.0227 


1.1884 


1.2003 


17 - 0000 


0.9360 


0.8618 


0.4918 


17 - 1200 


0.5325 


0.5421 


0.5940 


18 - 0000 


0.5678 


0.5863 


0.5923 


18 - 1200 


0.6549 


0 . 6600 


0.5514 


19 - 0000 


1.1639 


1.1672 


1.6403 


19 - 1200 


2.2758 


2.2377 


1.9609 


20 - 0000 


1.0744 


0.9818 


0.6177 


20 - 1200 


0.5083 


0.5216 


0.5058 


21 - 0000 


0.4897 


0.4766 


0.5948 


21 - 1200 


0.8447 


0.8370 


0.9129 


22 - 0000 


0.9807 


1.1674 


2.4002 


22 - 1200 


1.2914 


1.1747 


1.1335 


23 - 0000 


1.5845 


1.3996 


0.9415 


23 - 1200 


3.0389 


2.8514 


2.1848 
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TABLE XX 



Table of precipitable water vapor values (W Qj W cs , W cm ) 
for San Diego, Calif. 



DAY 

DAY - GMT 


OBSERVED 

w 0 


PROFILE 

W C s 


MULT. REG. 
Wcm 


16 - 0000 


0.7488 


0.7849 


0.5459 


16 - 1200 


0.8116 


0.8317 


0.6233 


17 - 0000 


0.8222 


0.8833 


0.6292 


17 - 1200 


0.8586 


0.8796 


0.8950 


18 - 0000 


0.8783 


0.8523 


0.7305 


18 - 1200 


1.0573 


1.0334 


0.7619 


19 - 0000 


0.9586 


1.0137 


0.8904 


19 - 1200 


1.1637 


1.2238 


1.6504 


20 - 0000 


0.7991 


0.9855 


1.5071 


20 - 1200 


0.7310 


0.8049 


0.8771 


21 - 0000 


0.6700 


0.7594 


0.8141 


21 - 1200 


0.9416 


0.9233 


0.8065 


22 - 0000 


0.7205 


0.8268 


1.0792 


22 - 1200 


1.2398 


1.3404 


1.8782 


23 - 0000 


1.1738 


1.3048 


1.7369 


23 - 1200 


1.7495 


1.7309 


1.4738 
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table XXI 



Table of precipitable water vapor values (W G , W cs , W cm ) 

for Spokane, Wash. 



DAY 

DAY - GMT 


OBSERVED 

W 0 


PROFILE 

w cs 


MULT. REG. 
Wcm 


16 - 0000 


0.3970 


0.3932 


0.2822 


16 - 1200 


0.3309 


0.3313 


0.2976 


17 - 0000 


0.3965 


0.4002 


0.4135 


17 - 1200 


0.3807 


0.3731 


0.6097 


18 - 0000 


0.3491 


0.3174 


0.2582 


18 - 1200 


0.4091 


0.3976 


0.3855 


19 - 0000 


0.4692 


0.4820 


0.4334 


19 - 1200 


0.5359 


0.5436 


0.5835 


20 - 0000 


0.5297 


0.5337 


0.6059 


20 - 1200 


0.6018 


0.4964 


0.3805 


21 - 0000 


0.4113 


0.3489 


0.3469 


21 - 1200 


0.3094 


0.2940 


0.1853 


22 - 0000 


0.3098 


0.3036 


0.4224 


22 - 1200 


0.7075 


0.7759 


0.9983 


22 - 0000 


0.7988 


0.7595 


0.6473 


23 - 1200 


0.9039 


1.0635 


1.3750 
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TABLE XXI 



Table of 


precipitable water 
for Sault Ste 


vapor values 
Marie, Mich. 


(Wo, w cs , W cm ) 


DAY 


OBSERVED 


PROFILE 


MULT. REG. 


DAY - GMT 


W 0 


W cs 


Wcm 


16 - 0000 


1.1656 


1.1476 


1.0986 


16 - 1200 


0.3112 


0.3207 


0.3864 


17 - 0000 


0.3320 


0.3293 


0.1813 


17 - 1200 


0.3511 


0.3401 


0.3188 


18 - 0000 


0.3865 


0.3693 


0.3514 


18 - 1200 


0.4089 


0.3887 


0.2807 


19 - 0000 


0.5890 


0.5605 


0.7646 


19 - 1200 


0.7004 


0.7405 


0.8876 


20 - 0000 


0.7978 


0.7349 


0.6034 


20 - 1200 


0.5112 


0.5162 


0.4981 


21 - 0000 


0.5046 


0.5150 


0.4953 


21 -- 1200 


0.2910 


0.2889 


0.4080 


22 - 0000 


0.3454 


0.3433 


0 . 4346 


22 - 1200 


0.2752 


0.2701 


0.2942 


23 - 0000 


0.3232 


0.3234 


0.3215 


23 - 1200 


0.2074 


0.2032 


0.1921 



59 



- Wes 
■* Wem 




60 



Waycross Ga . 



- Wes 




<N 

( UJ:> ) H31VM 3iaVildD3dd 



61 



Figure 4: Plot of W Q values, W cs values, and W cm values over the period 16 - 23 March, 1971 

for San Diego, Calif. 
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Figure 5: Plot of W Q values, w cs values, and W cm values over 

for Spokane, Wash. 
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Figure 6: Plot of W c values, W cs values, and W cm values over the period 16 - 23 March, 1971 

for Sault Ste. Marie. 



regression procedure for the date and station considered was A = 2.6793. 
Since the latter value is considerably larger than the sounding best-fit 
value of A , the application of the multiple regression A -value in 
Eq. (13) gives a much larger value of W than W ; hence the large error, 
W cm ~ W Q , at this sounding time (Figure 3). 

The results of Tables XIX, . . . ,XXII (and the four supporting diagrams 
Figures 3,..., 6) have been presented in the statistically summarized form 
as Table XXIII. In Table XXIII, the mean W Q -values have been compared 
to the computed mean values of W cg and W using the 16 sounding samples 
at the four stations, Waycross, San Diego, Spokane, Sault Ste Marie. 

The results were by-products of the correlation of W Q separately against 
W cs , and of W Q with W cm , which was applied to the two sets of data in 
Tables XIX, . . . ,XXII , respectively. 

The standard deviations of W Q , W cs , and W cm have been summarized in 
Table XXIII (column 3). Note that the means and standard deviations of 
the W-samples computed by the different procedures are nearly equal. 

The correlation coefficients R(W 0 , W cS ) and R(W Q , W cm ) are listed for 

the four stations in column 4 of Table XXIII. The corresponding 

2 

coefficients of determination R are listed in column 5, and are to be 

interpreted as the fractional explained variance (of W Q ) . Finally the 

column ”Std, Error” is to be interpreted as the mean dispersion of the 

curve of W Q minus W cs , and also that of W Q minus W cm , using the following 

✓ 

equation for standard error: 

Sw o |W c = (1 - R*)% a W 0 (4-8) 

Note that the estimation of W cs to give W Q leads to a standard error (a 
mean residual) of 0.07 cm, whereas the corresponding standard error in the 
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Statistics comparing the observed precipitable water values 
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Column Means .9680 .6495 .0701 .3093 



comparison of W cS with W Q is close to 0.31 cms. This larger error arises 
primarily from the southerly stations in Table XXIII. 

In conclusion, the value of W cs based upon a best-fit X -value in 
conjunction with Eq. (13) gives a good estimation of W^. On the other hand, 
the computation of W cm wherein the X -value is estimated by a multiple 
regression formula using gross parameters from the sounding has some skill 
in estimation of W Q , but perhaps needs additional study for selection of 
input variables for optimal inclusion in the multiple regression equation 
used to determine X . 



66 



VI. CONCLUSIONS 



At a group of 12 geographically diverse stations, the mixing-ratio 
profile was parameterized as a power-law profile, Eq. (6), having the 
characteristic exponent A . For each sounding A was derived by a 
least-squares fit with an accuracy usually characterized by a correlation 
coefficient of 0.95 or higher. At each station, the fit was most accurate 
in this form when the standard deviation of observed mixing-ratio was 
large. The set of A -values were grouped together in time series at each 
station, and they were found to be quite variable in general. This 
variability of A was related to the detailed structure of the moisture 
profile from the sounding data. A sizeable part of the standard devia- 
tion o' ( A ) was found to be explainable by use of a multiple regression 
procedure which incorporated easily defined properties of the sounding. 

The primary test of the A -profile parameterization was made in application 
to the comparison of observed versus computed values of the precipitable 
water vapor at four stations, with the latter given by 



of the observed precipitable water vapor value. Over a period of time 



gave moderately good verification when applied for the appropriate sound- 
ing using Eq. (13). There is some possibility that this procedure can be 
employed in a predictive sense provided the parameters used in the 
regression can be successfully generated by a primitive equations model. 




(13) 



On an individual sounding basis the A -value gave good verification 



it was possible to generate multiple regression equations for A which 
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It is possible that the multiple regression sounding parameters may be 
selected to yield closer agreement with the X -values determined from 
the individual soundings. 
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